arXiv: 1509.00975v2 [astro-ph.EP] 26Jun2016 


New members of Datura family 


Rosaev A.*^, Plavalova E.'^ 

°‘NPC Nedra, Yaroslavl, Russia 
^Astronomical Institute, Slovak Academy of Science, Slovakia 


Abstract 

Asteroid families are groups of minor planets which have a common origin in catastrophic disruptions. Young asteroid 
families are very interesting because they represent the product of their parent body’s fragmentation before orbital 
and physical evolutionary processes could have changed them. A group of minor asteroids associated with the largest 
body Datura(1270) is of particular interest because it has enough known members and resides in the inner part of the 
main asteroid belt and is easy to observe. Up to now, 7 members of this family are known. Here we discuss three 
new members of the Datura Family: 338309(2002 VR17), 2002 RH291 and 2014 OE206. To prove that these recently- 
discovered members belong to the Datura family, we conducted numeric integration with all gravitational perturbation 
over the last 800 kyrs. In the results, we have found that 338309(2002 VR17) and 2002 RH291 are very close to the 
mean orbit of this family throughout the calculation. In the case of 2014 OE206, it has a strongly chaotic orbit. The 
possible explanation of this is in the resonance character of its orbit. 
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1. Introduction 

The observations of the main asteroid belt have brought 
the discoveries of asteroid families groups of fragments 
from collisions of the large asteroids. Up to now, more 


than 120 asteroid families have been discovered |Nesvorny| 
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2006). The younger an asteroid family is, the more 


straightforward it is to analyse because they represent the 
product of their parent body fragmentation before orbital 
and physical evolutionary processes have had change them. 


2. Datura family 


A group of minor planets associated with the largest 
body Datura(1270), is of particular interest because it has 
enough known members and resides in the inner part of the 
main asteroid belt, which makes observation clearer. The 


Datura family was discovered by Nesvorny et al. (20061, 
after which, both its dynamical evolution and physical 
properties were rigorously studied. The model of rota¬ 


tional fission systems is described by (Pravec et al., 20101. 


According to their results, if the mass ratio of fragments 
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is significant, the rotational period of the primary grows 
longer. The very long period of rotation and the very elon¬ 


gated shape in the case of 2003 CL5(90265) (Vokrouhlicky 


et al. 20091, may be indirect evidence of its recent breakup 
and may highlight its possible mechanism. The Datura 
family consist of one large 10 km-sized asteroid (parent 
body) and a few small minor planets; possible fragments 
of a catastrophic breakup. Up to now only 7 members 


of the Datura family were known (Nesvorny et al. 2015 


2006 Vokrouhlicky et al.[ 2009), all of them belonging to 


the so-called S-type of minor planets Tholen (1984). 


3. The three new members of the Datura family 


As it is noted by Nesvorny et al. (2006), the discovery 


of new members of the Datura family is of great interest in 
order to study their group dynamics. Here we report three 
new members of the Datura family. In order to find new 
members of the Datura family, we studied 634109 orbits of 


asteroids from the Lowell Observatory catalogue (Howell 


et al. 1994). We extracted the bodies which possess the 


values of the orbital elements in intervals of 2.234 < a < 
2.236, 0.2 < e < 0.21, 5.9 < i < 6.1 according to the values 
of the orbital elements of 1270 Datura. We tested these 
bodies and several of their clones by numerical integration. 

In our research, we have found three new members of 
the Datura family in addition to those listed in |Nesvorny| 
et al. (2015 2006); Vokrouhlicky et al. (2009). These three 


new members are SSSSOff 2002 VR17, 2 002 RH291 and 
2014 OE206 (Rosaev & Plavalova, 2015). Tableshows 
the osculating elements for both the seven old members 
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Table 1: Osculating orbital elements of the Datura family members at epoch 16-01-2009 (JD 2454848). 


Number 

Object 

UJ 

n 

% 

e 

a 

1270 

Datura 

258.84072020 

97.88191931 

5.98954802 

0.207917356 

2.234316251 

60151 

1999 UZ6 

260.57151084 

96.79987218 

5.99341398 

0.207812911 

2.235186383 

89309 

2001 VN36 

266.85629891 

92.976202113 

6.02164403 

0.20634112 

2.23559504 

90265 

2003 CL5 

261.84375257 

95.69777810 

5.99543662 

0.207471114 

2.234865786 

203370 

2001 WY35 

260.44606900 

96.87210762 

5.99154981 

0.207435802 

2.235226671 

215619 

2003 SQ168 

259.39568050 

97.46790499 

5.99082425 

0.207908027 

2.234266272 

338309 

2002 VR17 

260.60046232 

96.80636247 

5.98940220 

0.207747446 

2.235169351 


2003 UD112 

263.12446374 

95.47877895 

6.00471739 

0.206934529 

2.234566113 


2002 RH291 

262.04302362 

95.75928088 

5.99579553 

0.207613218 

2.235326612 


2014 OE206 

261.72701226 

96.26530222 

6.00041657 

0.206975298 

2.235606600 


NOTE: The three new members of the Datura family are shown in bold. 


and the three new shown in bold. These new members 
of the Datura family are faint and small: as expected by 
absolute magnitudes, they are from 80 (2014 OE206) up to 
500 meters (2002 RH291 and 338309) in diameter, making 
them comparable with known member 203370 (Tab. [^. 
It is easy to see, that the orbits of the two new members 
338309(2002 VR17) and 2002 RH291 are determined much 
better than the orbit of 2003UD112, which was initially 


included in the Datura family (Nesvorny et al. 2006). And 


the precision of the orbit of 2014 OE206 is comparable to 
2003UD112 (Tab.[2|. 

To prove that these new members belong to the Da¬ 
tura family, we studied their numerical integration. The 
equations of the motion of the systems were numerically 
integrated 800 kyrs into the past, using the N-body inte¬ 


grator Mercury ([Chambers 1999) and the Everhart inte¬ 


gration method ([Everhart 19851. The initial state vectors 


for planets were obtained from JPL DE405 ephemerides 
at the same time as the elements of the Datura-group as¬ 


teroids (Standish et al. 1995). In addition, we use very 


few clones to establish the initial element’s uncertainties 
for each of the Datura familys new member’s orbital evo¬ 
lutions. The distribution of Datura family members in 
the coordinates (a, e) and and (z,e|^is given in Figure 
where the new members are shown as red dots and the 
original members as blue dots. All the asteroids of the 
Datura family, including new members, form a compact 
group in the orbital element’s phase space, true at least 
for 2002 VR17 and 2002 RH291 throughout all considered 
time intervals. In Figure [^ the evolution of the angu¬ 
lar elements for the two new members of the family are 
given, relative to 1270 Datura’s angular elements. It is 
evident that the values of the angular elements of the new 
members are close to the largest member of the family in 
all considered time intervals. For other orbital elements 
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Figure 1: The positions of the Datura family members in a, e (top 
figure) and i,e (botom figure) coordinates. The new members are 
shown as red dots and the original members as blue dots. 


'^Throughout the article we have used standard notations for or¬ 
bital elements a — semi-major axis in au, e — eccentricity, i - inclina¬ 
tion, Q - longitude of ascending node, ui — perihelion argument, tt — 
longitude of perihelion (the angular elements are in degrees). 
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Table 2: Some observational parameters of the new members of the Datura family in comparison with the originally known member 
2003 UD112. 


Number 

Object 

H 

Arc [days] 

1270 

Datura 

12.5 

37118.1 

60151 

1999 UZ6 

16.1 

6308.05 

89309 

2001 VN36 

16.2 

7915.94 

90265 

2003 CL5 

15.9 

7245.02 

203370 

2001 WY35 

17.1 

6828.77 

215619 

2003 SQ168 

17.0 

4702.05 

338309 

2002 VR17 

17.7 

5794.84 


2003UD112 

17.9 

25 


2002 RH291 

17.9 

4743.76 


2014 OE206 

19.3 

26 


NOTE: Some observational parameters of the new members of the Dat ura family. Where; H is the absolute magnitude, Arc is the orbital 
arc in days, spanned by observations used in an orbital computation by |Bowell et al.| l |1994[ |. 


we have similar relations. The orbital elements also re¬ 
main similar throughout the whole interval of integration 
(Fig. [31). The age of the Datura family is estimated to be 


that is to say that a question still lingers over the most 
precise age estimate. This inconsistency could mean that 
the question of the Datura family’s age is still open. 

As it seems in Figurej^ the orbital evolution of the new 
members of the Datura family is not in conflict with the 
estimated age of [Vokrouhlicky et"^ (2009) being 450 ± 
50 kyrs. The argument of the perihelion is somewhat bet¬ 
ter than the node longitude for the determination of age. 
In the case of 2014 OE206, we have a strongly chaotic 
orbit. According to our results, the semi-major axis of 
2014 OE206 significantly increased about 50 kyrs ago (Fig. 
[^. This effect presents a small variation of the initial data 
of this asteroid in the range of 2.23485 < a < 2.2360 au 
for the semi-major axis and at least 36.79° < M < 36.81° 
for the range of the mean anomaly. In fact, the orbit of 


450 ± 50 kyrs old according to 

Vokrouhlicky et al. 

(2009 

) 2.237 

and 530 ± 20 kyrs old according to ( 

Nesvorny et al. 

2006 

; 
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Figure 3: The positions of the Datura family members in a and e 
coordinates 300 kyr ago. Red point - 2014 OE206. 
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Figure 2: The average of differences in the mean longitude of 
the ascending nodes for 2002 RH291 (brown line) and for 
2002 VR17(338309) (green line); both belonging to 1270 Datura. 
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Table 3: Proper orbital elements of the original and new members of Datura family and the averaged over 300 kyrs integration semi-major 
axis(a) 


Number 

Object 

Average a [au] 

Synthetic proper elements 
a [au] e sin(f) 

1270 

Datura 

2.234694 

2.23469 

0.153413 

0.092307 

60151 

1999 UZ6 

2.234867 

2.23489 

0.153397 

0.092331 

89309 

2001 VN36 

2.235201 

2.23584 

0.153883 

0.092586 

90265 

2003 CL5 

2.234943 

2.23494 

0.153505 

0.092397 

203370 

2001 WY35 

2.234865 

2.23488 

0.153482 

0.092397 

215619 

2003 SQ168 

2.234663 

2.23465 

0.153450 

0.092318 

338309 

2002 VR17 

2.234959 

2.23494 

0.153503 

0.092313 


2002 RH291 

2.235002 

2.23490 

0.153533 

0.092316 


2014 OE206 is a typical example of so-called stable chaotic 
orbit (Milani et al., 19971. 

The possible explanation of this is in the resonance 
character of its orbit. As it was adopted by |Nesvorn^ 


et al. (2006), the orbit of 2001 VN36(89309) was assumed 


to be too uncertain because of resonance-related chaotic- 
ity (mean motion resonance 9 : 16 to Mars). Except 
noted above, there is a commensurability with a resonance 
of 7 : 2 to Jupiter in the vicinity of the Datura family, 
but its effect is weaker. The initial mean anomalies of 

2001 VN36(89309) and 2014 OE206 are very close, so we 
can suppose a finite range in the mean anomaly of the in¬ 
stability of the orbit caused by resonance. Resonance is 
also presented in the orbital evolution of the other mem¬ 
bers of the family but in a weaker form. The variations of 
the initial values of other orbital elements (eccentricity, in¬ 
clination, pericenter and node longitude) have a relatively 
small effect on the orbital evolution of 2014 OE206. 

Elementary consideration has shown that there has 
been no encounter with Mars closer than 0.01 au. How¬ 
ever, the change of semi-major axis in the case of 2014 
OE206 is significant; 50 kyrs ago, this asteroid leaped from 
one side of 9 : 16 resonance to another. On the other hand, 
its other orbital elements remained very close to the mean 
elements of the other members of the Datura family. So, 
we can conclude that 2014 OE206 has belonged to the Da¬ 
tura family for at least the last 50 kyrs. 

An additional argument which favours these new mem¬ 
bers belonging to the Datura family, is the calculation 
of the average orbital elements over the last 800 kyrs. 
Through numerical integration, we made a linear approx¬ 
imation of node longitudes. For 2002 VR17 (338309) and 

2002 RH291, we have a very close comparison in all the 
orbital elements. Proper elements, calculated by|Knezevic| 


& Milani (2003) confirm that 2002 VR17 (338309) and 
2002 RH291 belong to the Datura family (Tab. |^. In 
particular, we can point out a very close match of our cal¬ 
culations of averaged semi-major axis with proper values. 

A trait common in the Datura family is a linear de¬ 
pendence between Q and w (Fig. |^. This dependence 
is unique for the Datura family and remains unfound in 
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Figure 4: The positions of the Datura family members in coor¬ 
dinates. The new members are shown as red dots and the original 
members as blue dots. 


other young asteroid families (e.g. 1992YC2, Emilkowal- 
ski, etc Nesvorny & Vokrouhlicky (2006)), and studying 


it may be key to successfully reconstructing the dynamic 
history of the Datura family. 


4. Encounters in the Datura family 

As noted above, if we accept the hypothesis about the 
origin of all the members of the Datura family in the same 
event, the new members of the Datura family’s orbital 
evolution does not conflict with the sugested age of 450±50 
kyrs. 

Here we give a few arguments in favour of the sub¬ 
sequent breakup. Numerical integration shows, that two 
new members of the Datura family have a few recent close 
encounters with other known members. However, multi¬ 
ple encounters are far more interesting. Below we consider 
one interesting case of such an event. 

A true anomaly of asteroids is most sensitive at the 
semi-major axis of initial variations. If we speak about a 
difference of 10“® au in the initial semi-major axis, we have 
a variation in the asteroid position of about 10—15 degrees 
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Table 4: Distances in au between members of the subgroup close to epoch T = 157735 years ago. 


T,yr/d,au 1270 Datura / 1270 Datura / 

2003 SQ168(215619) 2001 VN36(89309) 

157730 0.0115 0.0160 

157740 0.0105 0.0207 


215619 / 

2001 VN36(89309) 
0.0178 
0.0145 


after 100 kyrs. In reality, this means that we cannot deter¬ 
mine the position of any minor planet along its orbit with 
the required degree of accuracy after a few thousand years. 
After such a time interval, we cannot determine the exact 
longitude of the studied close encounter and minimal dis¬ 
tance between the two asteroids. The YORP effect carries 
a strong influence in the evolution of semi-major axis; the 
maximum estimated drift rates of kilometer-sized objects 
being, da/dt = 2* . The mean anomaly can 

rotate with respect to an orbit with a hxed semi-major axis 
by 360° in 200 kyrs (Nesvorny & Vokrouhlicky 20061. 


Another point is, for young compact asteroid families, 
originating in low-velocity breakups, it is true that two 
minor planets spent a long time close to one another before 
impact. So, if we want to determine the epoch of a young 
asteroid family’s origin, we need to look for other criteria. 

As noted above, multiple encounters are far more inter¬ 
esting. Detection of such an event would be direct evidence 
of a breakup epoch. Here we note a triple close encounter 
with the participation of 1270 Datura, 2001 VN36(89309) 
and 2003 SQ168(215619) (Tab. (Hereinafter referred to 
as the subgroup). According to data from Tab. there 
is a small possibility that asteroid 2003 SQ168 (215619) 
originated in more recent event (Vokrouhlicky et al., 20091. 


Below we give some new suggestions in favour of this point 
of view. 

The data in Tab. |4] shows that three members of the 
Datura family were in very close proximity, with a distance 
smaller than 0.005 au of their distance from the Sun at 
an epoch of about 157735 ± 5 years from now. We can 
illustrate this fact as follows: 


T = 


Tio -f Ti 


13 


■Ts 


23 


6 Tij = Tij - T 


( 1 ) 


Denote as an epoch of the minimal distance between 
i and j asteroids, T as the mean epoch of the triple en¬ 
counter and 5Tij as the time interval at each partial en¬ 
counter to T. The values of 6 Tij at epoch 157735 ±5 years 
ago are by order of magnitude, smaller than other triple 
close encounters for these asteroids (Fig. [^. 

The differences in the averaged semi-major axes are 
determined by averaging their mean motion ie. the pe¬ 
riod between subsequent close encounters with other minor 
planets. In our case, this difference is about 5 x 10“® au. 
As it seems in Tab. the mean value of such differences 
in the Datura family is far more than the initial errors in 
the orbital elements. 


1.4 
1.2 
1.0 
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Figure 5: The differences in the period concerning these three bod¬ 
ies is illustrated as follows; red points are for 1270 Datura and 
215619 (2003 SQ168/306), blue points for 1270 Datura and 89309 
(2001 VN36) and green points for 215619 (2003 SQ168) and 89309 
(2001 VN36). 
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The period between two subsequent close encounters 
between a pair of asteroids in the family can be estimated 
using: 


P = 


P 1 P 2 

IP1-P2I 


\a^3/2 _ a^3/2\ 


( 2 ) 


Calculations based on Tab. 2003 SQ168(2156I9) P = 
160.54 ± 10 kyrs is in close agreement with our numerical 
integrations (Fig.|^. Between 1270 Datura and 2001 VN36 
(89309) we have a period of about 10 kyrs (9822 ± 12 yrs) 
again, in close agreement with our numerical integration. 

So, we can conclude that the period between two subse¬ 
quent close encounters between two asteroids of the Datura 
family is stable and has relatively negligible initial orbit er¬ 
rors. Of course, we cannot say this of the true anomalies 
or longitudes of the encounter point; they are too sensi¬ 
tive to initial conditions, but we can say this of the epoch 
of triple encounters with an error margin of about ±20 
kyrs, which is comparable with the accuracy of present 
age Datura family estimations based on angular elements 
analysis. At least we can detect the time intervals when 
mutual close encounters of studied asteroids cannot take 
place. For example, we can completely exclude the time 
intervals 0-140 kyrs and 180-290 kyrs as a possible origin 
of 2003 SQ168 (215619). 

It is interesting that if our conclusion about a subse¬ 
quent breakup close to the epoch 157.7 kyrs is true, this 
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Table 5: Time-dependence of the constant of the subgroup in a 300 kyrs time interval. 


Number 

Object 

Average a [au] 

Uncertainties in a 

1270 

Datura 

2.234694 

l.lxl0-« 

89309 

2001 VN36 

2.235201 

7.7x10-® 

215619 

2003 SQ168 

2.234663 

2.2x10-’^ 


means that non-gravitational forces (Yarkovsky etc.) do 
not have a significant effect on minor planets in the con¬ 
sidered subgroup. In any case, it is incorrect to exclude 

2001 VN36(89309) from Datura family dynamical recon¬ 
structions and to reject the possibility of the more recent 
origin of 2003 SQ168(215619). 

5. Conclusions 

Three new members of the Datura asteroid family have 
been found. These new members can help to improve our 
estimation of the epoch and circumstances of origin. Two 
of these new members (2002 VR17(338309) and 2002 RH291), 
have relatively well-determined orbits and belong to the 
Datura family throughout a time interval of at least 800 kyrs. 
The third new member, 2014 OE206, has a chaotic reso¬ 
nance orbit like 2001 VN36(89309) and is characterised 
by a significant increase in semi-major axis about 50 kyrs 
ago. If we accept the hypothesis about a one-moment 
origin of the Datura family, the dynamical evolution of 

2002 VR17 (338309) and 2002 RH291 does not conflict 
with the previous estimation of the age of the Datura fam¬ 
ily (450±50 kyrs). But studying mutual encounters in this 
asteroid cluster gives some arguments in favour of a subse¬ 
quent breakup at the epoch about 157.7±0.1 kyrs ago. In 
this more recent event three minor planets 1270 Datura, 
2001 VN36(89309) and 2003 SQ168(215619) took part. 
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